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Abstract 

We discuss several mechanisms of heavy quark production in (real) 
photon-nucleon and (real) photon - (real) photon collisions. In par- 
ticular we focuse on application of the Saturation Model. We discuss 
how to generalize the formula from virtual photon - proton scatter- 
ing and analyze threshold effects. We discuss a possibility to measure 
the cross section for 77 — > 2c2c. In addition to the main dipole- 
dipole contribution included in a recent analysis, we propose how to 
calculate within the same formalism the hadronic single-resolved con- 
tribution to heavy quark production. At high energies this yields a 
sizeable correction of about 30-40 % for inclusive charm production 
and 15-20 % for bottom production. We consider a subasymptotic 
component to the dipole-dipole approach. We get a good description 
of recently measured 17(77 ~~ > ccX). Adding all possible contributions 
to e + e~ — > bbX together removes a huge deficit observed in earlier 
works but does not solve the problem totally. Whenever possible, 
we compare the present approach to the standard collinear one. We 
propose how to distinguish different mechanisms by measuring heavy 
quark-antiquark correlations. 
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1 Introduction 



The total cross section for virtual photon - proton scattering in the region of 
small x and intermediate Q 2 can be surprisingly well described by a simple 
parametrization [0J inspired by saturation effects related to nonlinear phe- 
nomena due to gluon recombination for instance. This is called Saturation 
Model (SM) in the literature. The very good agreement with experimental 
data can be extended even to the region of rather small Q 2 by adjusting 
an effective quark mass (m e //). The value of the parameter found from the 
fit to the photoproduction data is between the current quark mass and the 
constituent quark mass. At present there is no deep understanding of the 
fit value of the parameter as we do not understand in detail the confine- 
ment and the underlying nonperturbative effects related to large size QCD 
contributions. 

We shall try to extend the succesful SM parametrization to quasi-real 
photon scattering. In the present analysis we shall limit to the production of 
heavy quarks which is, in our opinion, simpler and more transparent for real 
photons. Here we can partially avoid the problem of the poor understanding 
of the effective light quark mass, i.e. the domain of large (transverse) size of 
the hadronic system emerging from the photon. 

It was shown recently that the simple SM description can be succesfully 
extended also to the photon-photon scattering ||. In order to better un- 
derstand the success of such a description even more processes should be 
analysed on the same footing in the same framework. 

The heavy quark production is interesting also in the context of a deficit 
of standard theoretical QCD predictions relative to the experimental data as 
observed recently for bottom quark production in proton-antiproton, photon- 
proton and photon-photon scattering. This must be better understood in the 
future because usually the charm production in photon-photon collisions is 
considered as a good tool to extract the gluon distribution in the photon (see 
for instance [Q]). 

In the present paper we discuss and analyze many details of the heavy 
quark production simultaneously in photon-nucleon and photon-photon scat- 
tering. In particular we quantify some new terms not included so far in the 
literature on this subject. We discuss the range of applicability of the SM 
parametrization for heavy quark production. We put emphasis on some open 
unresolved problems and propose ways of their resolving. 
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2 Heavy quark production in photon-nucleon 
scattering 

In the picture of dipole scattering the cross section for heavy quark-antiquark 
(QQ) photoproduction on the nucleon can be written in general as 

a lN _ QQ (W) = J d 2 pdz \^ Q (p,z)\ 2 cr dN (P,^W) , (1) 

where <$>t is (transverse) quark-antiquark photon wave function (see for in- 
stance j|, H) and adN is the dipole-nucleon total cross section. Because of 
higher-order perturbative effects as well as nonperturbative effects the latter 
cannot be calculated in a simple way. In the following, inspired by its phe- 
nomenological success we shall use the Saturation Model parametrization 
for (JdN- It phenomeno logically incorporates color transparency at r <C 1/Q S 
and saturation at r > l/Q s - Here Q s is the saturation scale. Because for 
real photoproduction the Bjorken-x is not defined we are forced to replace 
x by a new appropriate variable. The most logical choice would be to use 
gluon longitudinal momentum fraction x g instead of Bjorken x. This would 
lead to only minor modification of parameters to describe experimental data 
for virtual photon - nucleon scattering. In the following we take 



a dN (p,z, W) = (T 

where 



exp 



1 ( x Q 



P 2 



-A/2 



(2) 



Qo \XoJ 

Above we have introduced the quantity xq = Q< ^^""" , where Mqq(z) = 

jtjzrpj- The latter definition (prescription) will become clear after the discus- 
sion below. 

In Fig.l we show predictions of the SM for charm photoproduction. The 
dotted line represents calculations based on Eq.(|l|). The result of this cal- 
culation exceeds considerably the fixed target experimental data 0, ^ ^ 0, 
lT| , [TjJ. One should remember, however, that the simple formula ([!]) applies 



at high energies only. At lower energies one should include effects due to 
kinematical threshold. In the momentum representation this can be done by 
requiring: Mqq < W, where Mqq is invariant mass of the final QQ system. 
This cannot be done strictly in the mixed representation used to formulate 
the Saturation Model as here the information about heavy quark/antiquark 
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momenta is not available. An approximate way to include the obvious lim- 
itation for heavy quark production is to neglect momenta and include only 
heavy quark masses in calculating Mqq. This upper limit still exceeds the 
low energy experimental data. There are phase space limitations in the re- 
gion x g — > 1 which has been neglected so far. Those can be estimated using 
naive counting rules. Counting reaction spectators we get for our process an 
extra correction factor 

s sup = (i - x g y . (4) 

In our case of mixed representation we are forced to use rather xq instead 
of x g . Such a procedure leads to a reasonable agreement with the fixed 
target experimental data as can be seen by comparing the solid line and the 
experimental data points. 

The deviation of the solid line from the dotted line gives an idea of the 
range of the safe applicability of the Saturation Model for the production 
of the charm quarks/ ant iquarks. The cross section for W > 20 GeV is 
practically independent of the approximate treatment of the threshold effects. 
The only data in this region come from HERA. Here the Saturation Model 



seems to slightly underestimate the HI collaboration data [[13 . 

For comparison in Fig. 1 we show the result of similar calculations in 
the collinear approach (thick dash-dotted line). In this approach the cross 
section reads as 

(7 7N _qq(W) = J dx N g N {x N ^ 2 F ) a^ QQ (W) , (5) 
where W is energy in the jg system. The gluon distribution in the nucleon 



is taken from Ref. ||16||. In the calculation shown in Fig.|] we have taken 
y?p — = mq + p\ with p± being the heavy quark/antiquark transverse 
momentum. The traditional collinear approach gives steeper energy depen- 
dence in the energy region considered than the SM predictions. In order to 
test both approaches more, experimental data for different energies W > 20 
GeV are needed. 

The calculation above is not complete. For real photons a soft vector 
dominance contribution due to photon fluctuation into vector mesons should 
be included on the top of the dipole contribution considered up to now. The 
VDM component seems unavoidable in order to understand the behaviour 
of F% — F£ at small photon virtualities |T8| . Furthermore it seems rather 
a natural explanation of the strong virtuality dependence of the structure 
function of the virtual photon as measured recently pil . 
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In the present calculation we include only dominant gluon-gluon fusion 
component. In this approximation 



^n%q( w ) = % I dx vdx N gv(x v ,/i 2 F ) g N {xN,^ 2 F ) (J gg _ Q Q{W) . 

y JV J 

(6) 

Here the fy constants describe the transition of the photon into vector mesons 
(p, u, <fi) and are taken in the on-shell approximation from the decay of vector 



mesons into dilepton pairs []18f taking into account finite width corrections. 



The gluon distributions in vector mesons are taken as that for the pion [15| 



9v(xv, /4) = 9v( x *> /4) ( 7 ) 



and that in the nucleon from Ref.|]16[. For factorization scale we take fi 2 F = 
m Q + P 2 ± w ith Pi. being the heavy quark/ ant iquark transverse momentum. 
The dash-dotted line in Fig.^ shows the VDM contribution calculated in 
the leading order (LO) approximation for cr gg ^QQ (see e.g. JH]]). The so- 
calculated VDM contribution cannot be neglected at high energies. Taking 
into account that usually the next-to-leading order approximation leads to 
an enhancement by a factor K ~ 2, the VDM contribution is as important 
as the continuum calculated above. 

The situation for bottom photoproduction seems similar. In Fig.[| we 
compare the Saturation Model predictions with the data from the HI collab- 



oration HT3 |. Here the threshold effects may survive up to very high energy 
W ~ 50 GeV. Again the predictions of the Saturation Model are slightly 
below the HI experimental data point. The relative magnitude of the VDM 
component is similar as for the charm production. 

The Saturation Model slightly underestimates high-energy data for both 
charm and bottom production. The parameters of the Saturation Model: 
cr and 77i e // are to some extend correlated when extracted from the fit to 
the experimental data. In principle, because in the case of heavy quark 
production one is free of the uncertainty of the effective quark mass, one could 
allow to modify <jq to describe the HERA data on heavy quark production. 
It is obvious that then one would need to in 
shall leave this option for a separate study. 



It is obvious that then one would need to increase m e // to describe crffl. We 



Heavy quark production in photon-photon 
scattering 
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3.1 Heavy quark-ant iquark pair production 

In the dipole-dipole approach (see for instance ||) the total cross section for 
77 ~~ > QQ production can be expressed as 

o%-> Q q(W) = Yl /l^(Pi^i)ri^ /2/ HP2,^)|V dd (p 1 ,p 2 ,a;Q / )dV^irf 2 P2^2 
+ Yl J \^ flIl (puZi)\ 2 \^ QQ (p2,z 2 )\ 2 a dd (p 1 ,p 2 ,x fQ ) d 2 p 1 dz 1 d 2 p 2 dz 2 



h+Q' 



where are the quark-antiquark wave functions of the photon in the 
mixed representation and add is dipole-dipole cross section. While for the 
heavy quark-antiquark pair the photon wave function is well defined, for 
light quarks one usually takes the perturbatively calculated wave function 
with the quark/antiquark mass replaced by m e ff. This parameter provides 
a useful cut-off of large-size configurations in the photon wave function. 

There are two problems associated with direct use of (§). First of all, it 
is not completely clear how to generalize the dipole-dipole cross section from 
the dipole-nucleon cross section parametrized in Secondly, formula (||) 
is correct only at sufficiently high energy W 2m,Q. At lower energies one 
should worry about proximity of the kinematical threshold. 

In a very recent paper a new phenomenological parametrization for 
the azimuthal angle averaged dipole-dipole cross section has been proposed: 



ah/ \ o 

(?dd [Xab,Pl,P2) = ^0 



exp 



r eff 



4i2§(zo&) 



(9) 



Here 



1 / \ ~ A /2 



Ro(x ab ) = — — (10) 
Qo V x J 

and the parameter x a b which controls the energy dependence was defined as 

_ 4m 2 a + 4m 2 b m) 
Xab- w2 ■ {H) 

In order to take into account threshold effects for the produtcion of qqq'q' an 
extra phenomenological function has been introduced [[| 

Sthresh(Xab) = (1 — X a b) 5 (12) 



(8) 
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which is put to zero if x a b > 1. This factor strongly reduces the cross section 
at low energies. Different prescriptions for r e ff have been considered in ||, 

2 2 

with r 2 ^ = pi+pi being the best choice phenomenologically ||. 

The phenomenological threshold factor (0) does not guarantee automatic 
vanishing of the cross section exactly below the true kinematical threshold 
W = 2m a + 2mb- Therefore instead of the phenomenological factor we rather 
impose an extra kinematical constraint: Mjj + Mqq < W on the integration 
in @. The use of extra dynamical factor (1— x a b) 5 for heavy quark production 
in photon-photon collisions is in our opinion disputable. Instead, in the 
present paper we shall estimate the effect of damping of the cross section 
in the neighbourhood of threshold due to simple kinematical limitations on 
the final quark/ ant iquark transverse momenta. As will be discussed in the 
course of this paper this leads to a similar suppression at least numerically. 

Identically as for 'yN scattering in the mixed representation, where the 
transverse momenta are not given explicitly, the quark- ant iquark invariant 
mass is not well defined. We suggest therefore to use rather the effective (z- 
dependent!) invariant masses < M q q > with transverse momenta neglected 
and effective quark mass as used in many other mixed representation calcu- 
lations for jN scattering known in the literature. It is not completely clear 
how to generalize the energy dependence of a^N m photon-nucleon scattering 
to the energy dependence in a^d in photon-photon scattering. In the follow- 
ing we define the parameter which controls the Saturation Model energy 
dependence of the dipole-dipole cross section in a symmetric way as 

x Qf = x Qf {z 1 ,z 2 ) = C-{<M QQ >{z 1 )+<M ff ->(z 2 )) 2 /W 2 , 

x fQ = x fQ (z 1 ,z 2 ) = C-(<M ff ->(z 1 )+<M QQ >(z 2 )) 2 /W 2 (13) 

instead of ([11]). Here we have made explicit the dependence on Z\ and z 2 . 
In the following we shall use C = 1 and only in some cases compare to the 
results with C = 0.5. [] f\ In Fig.|3| we compare our prescription with those 
used in for the cc production (left panel) and for the bb production (right 
panel). In comparison to J| our prescription leads to a small reduction of the 
cross section far from the threshold and a significant enhancement close to 
the threshold. Some consequences of this will be discussed separately in the 
context of the observed excess of bb production in positron-electron collisions. 
For comparison we show in Fig.|^ also a result obtained in the two-gluon 

1 There is only small difference between the two prescriptions. 

2 By construction < xq/, x/q < 1. 
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exchange model. In this case 



Mpi,p 2 ) = \J " (2 - - e-^) • (2 - - e**) . 

(14) 

In principle one can allow for running of a s (fi 2 ). However, the choice of the 
scale fi 2 is not completely clear. In the present calculation a rather large 
constant value a s = 0.35 was taken. If a s (m 2 /m 2 ) was used the cross section 
would be negligibly small. Thus it becomes clear that the Saturation Model 
leads to a huge enhancement relative to the two-gluon exchange model high 
above the threshold. Close to the threshold both results almost coincide. 
The departure of the 2g-exchange result from the constant value below W m 
3 • (2mq) is due to the threshold effects. 

In calculating the main SM component above we have considered only 
obvious kinematical limitations possible to implement in the mixed repre- 
sentation. As already discussed it is not possible to include the transverse 
momentum limitations in the mixed representation. Now we shall estimate 
the effect due to the kinematically limited integration over transverse mo- 
menta of final heavy quarks/ ant iquarks. This effect can be taken into account 
consistently only in the momentum representation. This would require a 
reformulation of the whole Saturation Model and clearly goes beyond the 
scope of the present analysis. In order to gain experience we have studied 
first the effect of the kinematical constraint on the transverse momenta in 
the two-gluon exchange approximation in the momentum representation. In 
this case we can easily both include and exclude the kinematical limitations 
on transverse momenta. Then one can approximately correct our mixed rep- 
resentation calculation as 

^QqWU = ^ d QQ {W)\ awroX ■ Rc/a{W) . (15) 

Here for brevity we have introduced the ratio: 



2g,mom (xx/\ I 
(T 11 -,QQ\ VV ) \approx 



R c /a(w) = 2 ;:„r;„„ . (ie) 



In Fig. | we show the so-obtained correction factor as a function of W for both 
cc (solid) and bb (dashed) production. As can be seen from the figure the 
correction is significant even far from the threshold, i.e. leads to a significant 
damping of the cross section. We shall discuss this damping in the context 
of the bb deficit in a separate section. 
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3.2 2Q2Q final states 

The dipole-dipole approach in general, and the Saturation Model as its par- 
ticular realization, leads to a unique prediction for the 2Q2Q (two identical 
heavy quarks and two identical heavy antiquarks) production in the final 

state 

°^2Q2q( W ) = j \^l Q {pl^l)\ 2 \^ Q {p2,Z 2 )\ 2 a dd { Pl ,p 2 ,X QQ ) d 2 Pl dz x d 2 p 2 dz 2 

(17) 

The same prescriptions are used here as for the QQ final states in the previous 
section. 

In Fig.|5] we compare our prescription to that from || for four heavy 
quark production for Q=c (left panel) and for Q=b (right panel). Here there 
is even larger difference between the two prescriptions than for the heavy 
quark-antiquark pair production. In Fig.|| we show the ratio defined as 

R^(W) = a ^ 2Q2Q ( Z ] ■ (18) 

Sufficiently above the kinematical threshold W ^> 4mQ the ratio saturates 
at the level of about 8 % and 1 % for 2c2c and 2626, respectively. The 
fluctuations of the order of 1 % that one can observe in the figure are of 
numerical origin. 

The predictions for 77 — > 2Q2Q shown in Fig.|5] and ||| are the only ones in 
the literature we know. In the standard collinear approaches the 2Q2Q final 
states can be produced only in next-to-leading order calculations or/and in 
the hadronization process, if the charmed mesons, e.g. D*, are measured to 
identify charm quarks/antiquarks. It would be interesting to compare in the 
future the present result with the standard collinear NLO predictions. Also 
from the experimental side the 2Q2Q production in photon-photon collisions 
is terra incognita. In our opinion, the 2Q2Q channel has a better chance to be 
a stringent test of the dipole-dipole approaches and the Saturation Model in 
particular. It is not clear to us if the experimental verification can be feasible 
with the present LEP2 statistics. It seems, however, possible with the future 
photon-photon colliders like that planned at TESLA (see for instance ||). 

3.3 Short- versus long-distance phenomena 

What are typical distances probed in heavy quark production ? Is the heavy 
quark production a short distance phenomenon ? These questions can be 
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easily answered in the mixed representation formulation considered in the 
present paper. In Fig.^j we display the integrand of 



The maxima in Fig.0 correspond to the most probable situations. For one 
light (m„ = ma = m , m s = m + 0.15 GeV) and second heavy quark- 
antiquark pair the map is clearly asymmetric. One can observe a ridge par- 
allel to the pi or p 2 axis. There is not well localized maximum. Both short 
and long distances are probed. 

For comparison, in the bottom part of the figure, we show similar maps 
when both pairs consist of light (u,d,s) quarks/antiquarks (left-bottom) and 
in the case when both pairs consist of charm quarks/antiquarks (right-bottom) 
For light quarks (u,d,s) one observes a clear maximum at (pi, p 2 ) = (1 GeV -1 , 
1 GeV -1 ) = (0.2 fm, 0.2 fm). In this case a nonnegligible strength extends, 
however, up to large distances p\ and p 2 . Only in the case of the production 
of two cc pairs the cross section is dominated exclusively by short-distance 
phenomena. 

3.4 Hadronic single-resolved processes 

Up to now we have calculated the contribution when photons fluctuate into 
quark-antiquark pairs. Then the final quark and antiquark carry away the 
whole longitudinal momentum of the parent photon and are predominantly 
emitted in the same photon hemisphere. In the standard collinear approach 
one usually includes so-called resolved contributions, when heavy quark- 
antiquark pairs are created either in the photon-gluon / gluon-photon fusion 
or in the gluon-gluon fusion and quark-antiquark annihilation processes. In 
the first case, known as the single-resolved process, only a small fraction of the 
first or the second photon longitudinal momentum enters into the production 
of heavy quark or antiquark. In the second case, known as double-resolved 
process, this is true for both photons. The arguments above demonstrate 
that the resolved processes are different, at least kinematically, from those 
included in the Saturation Model, or more generally in the dipole-dipole in- 
teraction picture. This means that the resolved processes are not included 
in the dipole-dipole scattering approach. At the same time the standard 
collinear approach is not complete too. Can the dipole approach be sup- 
plemented to include the hadronic resolved processes? Let us try to answer 
this question by combining a simple vector dominance model and the dipole 
picture. For simplicity we shall limit to the dominant photon fluctuations 
into light vector mesons: p, u and 0, only. It is sufficient to include the light 




(19) 
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vector mesons because the contribution of heavier mesons, as being highly 
off-mass-shell, is expected to be considerably suppressed. 

If the first photon fluctuates into the vector mesons the so-defined single- 
resolved contribution to the heavy quark-antiquark production can be calcu- 
lated analogically to the photon-nucleon case as 

= ^]fj \^ Q (P2,z 2 )\ 2 a Vld (p 2 ,x 1 ) d 2 p 2 dz 2 , (20) 

where $2 i s the second photon QQ wave function and ay 1( i is vector me- 
son - dipole total cross section. In the spirit of the Saturation Model, we 
shall parametrize the latter exactly as for the photon-nucleon case |]] with a 
simple rescaling of the normalization factor (Jq V = 2/3 ■ <Jq N . In the present 
calculation (Jq N as well as the other parameters of SM are taken from |L|. 
Fully analogically if the second photon fluctuates into vector mesons 

= EtJ / i^v^i)rw Pl ,* 2 ) d 2 Pl d Zl . (21) 

This clearly doubles the first contribution ( |20"D to the total cross section. 
Of course, if the rapidity distributions of heavy quark/antiquark are con- 
sidered, both contributions must be treated independently. We leave the 
analysis of (rj, k±) distributions for a separate study. The integrations in 
( |2T)| ) and fl2"T| ) are not free of kinematical constraints. When calculating both 
single-resolved contributions it should be checked additionally if the heavy 
quark-antiquark invariant mass Mqq is smaller than the total photon-photon 
energy W . As in the previous cases in the mixed representation it can be 
done only approximately. In analogy to scattering the most logical defi- 
nition of the parameter which controls energy dependence of <iy x d or <Jdv 2 i s: 
x 2 /i = Mqq(z 2 /i) /W 2 and corresponds to the gluon momentum fraction in 
the second or first photon (vector meson), respectively. 

In Fig|| we compare the present result with the result of the leading order 
collinear approximation. In the latter case 



SR,coll ST^ 47T f . 2 . 

a ^ Q Q ~ 2^ IT / 9vAXVi, Pf) a gj^QQ ax Vi 



+ E 



Vl JVi 

4:71 



f2 

Jv 2 



T I 9v 2 (. x v 2 , Hf) Vjg^QQ dx V2 . (22) 

v 2 



Also in this calculation the gluon distributions in the vector mesons are taken 
as that for the pion [jnj g v (x g ,fip) = g^Xg, fi 2 F ). The scale p? F = m 2 c + p 2 ± 



11 



is taken in practical calculation with p± being the transverse momentum of 
the final heavy quark. 

There is a difference in shape of the cross section as obtained in the 
SM (thin solid) and collinear (dashed) approaches. Part of the difference 
may come from the fact that up to now we have not included phase space 
limitations when X1/2 — > 1 which may be important at low energies. In 
this case a naive application of counting rules would lead to the following 
suppression factor 

Ssr = (1 - x 1/2 ) 2n ^- 1 . (23) 

Following standard prescription by making use of the picture that the un- 
resolved photon probes the resolved one and assuming that a vector meson 
consists of a valence quark and antiquark we get usr = 3. The SM result 
corrected by (^) is shown by the thick solid line in Fig.§. After the correc- 
tion the results obtained from the two different approaches are numerically 
fairly similar. 



3.5 Summing all contributions to inclusive cross sec- 
tion 

In Fig.y we show different contributions to the inclusive c/c (left panel) and 
b/b (right panel) production in photon-photon scattering. The thick solid 
line represents the sum of all contributions. 

Let us start from the discussion of the inclusive charm production. The 
experimental data of the L3 collaboration [2(J are shown for comparison. The 



modifications discussed above lead to a small damping of the cross section in 
comparison to ||. The corresponding result (long-dashed line) stays below 
the recent experimental data of the L3 collaboration [EO]. The considered 
in the present paper hadronic single-resolved contribution constitutes about 
30 - 40 % of the main Saturation Model contribution. At high energies the 
cross section for the 2c2c component is about 8 % of that for the single cc 
pair component (see Fig.|5]). In inclusive cross section its contribution should 
be doubled because each of the heavy quarks/antiquarks can be potentially 
identified experimentally. In principle the events with 2c2c can be subtracted 
both when flavour tagging is applied or charmed mesons are measured to 
identify c or c. In practice, because the efficiency of the flavour tagging is 
very small and only a small fraction of 2c2c can be removed, the measured 
inclusive cross section is two times bigger than the integrated cross section 
to the 2c2c final state. 
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At higher photon-photon energies the direct contribution is practically 
negligible. This is in contrast to the energy dependence in the positron- 
electron collisions. Here even at large energies the direct contribution consti- 
tutes almost half of the corresponding cross section. The reason is that even 
at high e + e~ energies the contributions of small photon-photon energies are 
dominant which can be easily understood in the equivalent photon approx- 
imation to be discussed in the next section. The hadronic double-resolved 
contribution, when each of the two photons fluctuates into a vector meson, 
calculated as 

a j^gQ = ^2 JTJT I 9Vi(xvi,Vf) 9v 2 0ev2>/40 a gg->Qq(W) dx Vl dx V2 

Vi V2 Vl V2 

(24) 

and shown by the thin solid line in the figure becomes important only at 
very high energies relevant for TESLA. Here we have consistently taken 

g v (x v , /4) = g*(x v , H%) and /4 = m Q + P±- 

The situation for bottom production (see right panel) is somewhat dif- 
ferent. Here the main SM component is dominant. Due to smaller charge 
of the bottom quark than that for the charm quark the direct component 
is effectively reduced with respect to the dominant SM component by the 
corresponding ratio of quark/antiquark charges: (1/9) 2 : (4/9) 2 = 1/16. The 
same is true for the 2626 component. Here in addition there are threshold 
effects which play a role up to relatively high energy. Also the single-resolved 
component is relatively smaller which has no simple explanation. 

The role of different mechanisms for charm and bottom production is also 
summarized numerically in Table 1, where we have collected corresponding 
cross sections for a few selected values of photon-photon energies. 



3.6 e + e- -» e+e"66X 

Up to now no attempt was done to unfold experimentally the cross section 
for 77 — > bbX. Only cross section for the e + e~ — > bbX reaction was obtained 



recently by the L3 and OPAL collaborations at LEP2 [ 20[ ^2| . The measured, 
both positron and electron antitagged, cross sections can not be described 
as a sum of direct and single-resolved contributions, even if next-to-leading 



order corrections are included f22fl . The measured cross section exceeds the 
theoretical predictions by a large factor. This is a new situation in comparison 
to the charm production where the deficit is much smaller. We shall consider 
the case of bottom production in more detail below. 

The cross section for the e + e~ — > bb X reaction when both positron 
and electron are antitagged can be easily estimated in the equivalent photon 
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approximation (EPA) as 



cr(e + e~ -> bbX; W ee ) = J dx A dx B /a (#6, ^a)/b(^6, 

(t( 7 7 -> W) , (25) 

where /a and are virtuality-integrated flux factors of photons in the 
positron and electron, respectively, and 9 max is the maximal angle of the 
positron/electron not to be identified by the experimental aparatus. In the 
present analysis we calculate the integrated flux factors f\ and fs in a sim- 
ple logarithmic approximation. The photon-photon energy can be calculated 
in terms of photon longitudinal momentum fractions xa and Xb in positron 
and electron, respectively, as W = ^XAXss e e- It is instructive to visualize 
how different regions of W J7 contribute to cr(e + e~ — ► bbX;W ee ). For this 
purpose it is useful to transform variables from xa, xb to xp = xa — %b and 
W 77 = W. Then 

a(e + e _ -> W ee ) = ^ Gfe F cW J" fA(E b , 9 max ,x A )fB(E b , 9 max , x B ) 

(7(77 -> WX; W) , (26) 

where the jacobian J is a simple function of kinematical variables. 

Before we go to the bb production let us look at the corresponding cc pro- 
duction. In Fig.|Il] we compare the x ir-integrated cross section da j dW(e + e~ — 
ccX) with recent experimental data of the L3 collaboration p0| . Different 
mechanisms are shown separately. The agreement with the L3 collaboration 
experimental data here is even better than for the unfolded 77 — > cc data 
shown in Fig.^. This is probably due to our simplified treatment of the pho- 
ton flux factors. The good agreement of the sum of all contributions with 
the data gives credit to our model calculations. 

The integrand I(xp,W) of (^6|) for bb is shown in Fig.|ll] for the direct 
production (left panel) and for the Saturation Model (right panel) including 
all contributions considered in the present analysis. Quite different pattern 
can be observed for both mechanisms. While for the direct production one 
is sensitive mainly to low-energy photon-photon collisions, in the Saturation 
Model the contributions of high-energies can not be neglected and one has 
to integrate over essentially up to W ee . This difference in W 77 is due to 
different energy dependence of (1(77 — > bb; W 77 ) for the different mechanisms 
considered as has been discussed above. But even in the latter case the 
integrated cross section is very sensitive to the region of not too high-energies 
W ~ 20 GeV, where the not fully understood threshold effects may play 
essential role. 
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For LEP2 averaged energy < W ee >~ 190 GeV the cross section inte- 
grated taking into account experimental cuts is a(e + e~ — > bbX) = 6.1 pb 
(C = 1) or 7.4 pb (C = 1/2) for the dipole-dipole SM scattering process. 
If the limitations on transverse momenta are included in addition through 
the factor R c / a (see (0)), then the corresponding cross section is reduced 
to 3.2 pb (C = 1) or 3.9 pb (C = 1/2). f] Corresponding cross section for 
the direct production is cr(e + e~ — > bbX) = 1.2 pb. The hadronic single- 
resolved contribution calculated here in the Saturation Model as described 
in the present paper is very similar in size to that calculated in the standard 
collinear approach |22]. As can be seen in Table 2 the 2626 contribution is 
practically negligible. We have completely omitted the double-resolved con- 
tribution which is practically negligible (see Table 1). The sum of the direct, 
66 SM, 2626 SM and the hadronic single resolved SM component is 9.3-10.6 
pb in the case when no transverse momenta cuts on the main SM component 
are included and 6.4-7.1 pb with the cuts. These numbers should be com- 
pared to experimentally measured cr(e + e~ — > bbX) — 13.1 ± 2.0 (stat) ± 2.4 
(syst) pb || (L3) and preliminary a{e + e- bbX) = 14.2 ± 2.5 (stat) ± 5.0 
(syst) pb [p3] (OPAL). In comparison to earlier calculations in the literature, 
the theoretical deficit is much smaller. The success of the present calculation 
relies on the inclusion of a few mechanisms neglected so far, in particular 
the dipole-dipole contribution which, in our opinion, is not contained in the 
standard collinear approach. 

Up to now only W-^-integrated cross section has been determined exper- 
imentally. This, in fact, does not allow to identify experimentally whether 
the problem is in low or high W 7J . In order to identify better the region 
where the standard collinear approach fails it would be useful to bin the ex- 
perimental cross section in the intervals of W 77 making use of a possibility 
to measure W V i s which can be related to W 77 via a suitable Monte Carlo 
program. At present, even splitting the cross section for e + e~ — *> bbX into 
cr(W 77 < W ) and a(W 71 > W ) for W ~ 20 GeV would be useful and should 
shed more light on the problem of the experimental excess of 66 relative to 
the "standard" QCD approach. 



3.7 Subasymptotic components to the dipole-dipole ap- 
proach to 77 — > QQ 

The Saturation Model by construction includes only dominant asymptotic 
contributions relevant at very high energies. As discussed above the problem 

3 For comparison in the two-gluon exchange approximation with a s — 0.35 the corre- 
sponding cross section is 10.9 pb. 
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of heavy quark production (especially bb) may be a bit more complicated. In 
the following we shall try to generalize the dipole-dipole approach to include 
dynamics which may be of relevance also close to threshold. As discussed 
above this can be the region of the deficit of theoretical predictions observed 
for the bottom production in photon-photon collisions. 

In the following as an example we shall try to estimate the cross section for 
the process when the quark associated with one photon annihilates with the 
same-flavour antiquark associated with the second photon. Then the heavy 
quark-antiquark pair can be produced via exchange of s-channel gluon. In 
the dipole-dipole approach this effect can be estimated(!) in terms of the 
familiar quark-antiquark photon wave functions as follows 

<T%+QQ = £> l/2 [ d Zl d 2 Pl dz 2 d 2 p 2 \^ h (z 1 , Pl )\ 2 \^(z 2 ,p 2 )\ 2 

fuh 

■ [^qq-^QQ^, P2) + Vq q ^QQ{s, Pi, &)] ■ (27) 

The Kronecker 5 reflects flavour conservation in QCD. In the following we 
shall include only light flavours (u,d,s) as "constituents" of both photons 
and use the same treatment (suppression) of large size physics as in all cases 
before. For the purpose of the estimate a sensible approximation is to neglect 
light quark/ antiquark transverse momenta and write 

<T-00 = Ylj VfhM) ?ffj&) <Tqr+QQ{8) d Zl dz 2 

9//i<>2) ^-qq(s) dz x dz 2 , (21 

f 

where a qq -_,QQ or a qq _,QQ are calculated in the Born approximation with 
collinear incoming light quark and antiquark [17]. The argument of a s (p 2 R ) 
is set to p? R = s in practical calculations. In Eq.fl2~8|) we have introduced for 
brevity 

qf(z)= [\9"(z,p)\ 2 <fp. (29) 



E 



Because the perturbative photon wave function is singular at p = 0, the 
expression above is formally divergent. However, there is a natural cut-off 
at small p's: p m in ~ 1-/Vs which makes the integral finite. This cut-off 
corresponds to the upper limit on transverse momenta in the momentum 
representation calculations. 

In Fig.|I2| we show the cross section for the subasymptotic contribution 
considered for both cc and bb production. The maximum of the cross sec- 
tion is concentrated in the close neighbourood of corresponding kinematical 
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thresholds at W ~ 10 GeV and W ~ 30 GeV for charm and bottom, respec- 
tively. This concentration of the strenght at threshold justifies the introduced 
name "subasymptotic" for the process considered. The main (asymptotic) 
component of SM is shown for comparison by the dashed line and the QPM 
contribution by the dotted line. The subasymptotic contribution considered 
shows similar energy dependence as the QPM component. [] At high ener- 
gies it is rather small in comparison to the leading asymptotic component. 
Only, at threshold it constitutes a nonnegligible fraction of the heavy quark 
production cross section. When integrating crp^^x with virtual photon flux 
factors over W 77 and ij? we obtain <7^+ b e -^ b i x ~ 0.1 pb at W ee = 190 GeV. 
This is rather a conservative estimate because in the present calculations we 
have used leading order formulas for a q q^QQ. This seems tiny in compari- 
son to the contributions considered before (see Table 2). We could increase 
slightly the cross section by a different choice of the renormalization scale [i 2 R 
in calculating a qq _,QQ. We expect that in reality <7 s + 6 _ , ^ should not exceed 
0.5 pb. 



3.8 Quark-antiquark correlations 

So far mainly integrated cross section for heavy quark/ ant iquark production 
was considered in the literature. Only in a few cases inclusive distributions in 
transverse momentum or rapidity (see e.g.^4|) were presented. No attempts 
were done so far to analyze the final state in more detail. In our opinion inves- 
tigating correlations between heavy quark - heavy antiquark could be much 
more conclusive to identify the production mechanisms than the integrated 
cross section or even a single variable distribution. 

In principle any correlation between two kinematical variables of the final 
quark and antiquark would be of interest. We suggest that the following final 
quark/antiquark momentum fractions: 

_ PQ ~ 
XQ ~ \PQ\ Ul1 ' 



xq = prn 7l , (30) 
IPC 



where 



K = p^r (31) 

Kil 



4 In fact the contribution considered can be viewed as higher order correction to leading 
order QPM. 
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would be very useful to separate the different mechanisms (approches) anal- 
ysed in the present paper. In the definition above j>q and pq are momenta 
of heavy quark and antiquark, respectively, and p 7l is the momentum of the 
first photon, all in the photon-photon center of mass frame. By definition -1 
< xq,xq < 1. Similar quantities are being used at present when analyzing 
e.g. jet production at HERA to separate out resolved and direct processes. 

In Fig.|l3| we present a sketch of naive expectations. Although a precise 
map requires detailed calculations for each mechanism separately, which is 
beyond the scope of the present analysis, it is obvious that the separation of 
different mechanisms here should be much better then for any inclusive spec- 
tra. In the case of dipole-dipole approach (the elongated ellipses) this would 
require to go to the momentum representation. The mixed representation 
used in the present paper is useful only for integrated cross sections. 

Experimentally, the suggested analysis would be difficult at LEP2 because 
of rather limited statistics. We hope that such an analysis will be possible at 
the photon-photon option at TESLA. At present, even localizing a few LEP2 
coincidence cc events in the diagram x c versus x 5 would be instructive. 

4 Conclusions 

There is no common consensus in the literature on detailed understanding of 
the dynamics of photon-nucleon and photon-photon collisions. In this article 
we have limited the discussion to the production of heavy quarks simultane- 
ously in photon-nucleon and photon-photon collisions at high energies. 

A special emphasis has been put on the application of the Saturation 
Model which turned out recently very succesfull in the description of ex- 
perimental data for DIS at small Bjorken x. We have suggested how to 
generalize the model to applications with real photons. The sizeable mass of 
charm or bottom quarks sets natural low energy limit on naive application of 
the Saturation Model. Here a careful treatment of the kinematical threshold 
is required. In the mixed representation used to formulate the Saturation 
Model the effect can be included only approximately. 

We have started the analysis from (real) photon-nucleon scattering, which 
is very close to the domain of the Saturation Model as formulated in P]. If 
the kinematical threshold corrections are included the SM gives numerically 
similar results as the standard collinear approach for both charm and bottom 
production. We have numerically estimated the vector dominance contribu- 
tion to the heavy quark production. 

The major part of the present analysis has been devoted to real photon - 
real photon collisions. For the first time in the literature we have estimated 
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the cross section for the production of 2c2c final state. Furthermore we have 
discussed how to include this component to the inclusive charm production 
as derived in the present experimental analyses. We have found that this 
component constitutes up to 10-15 % of the inclusive charm production at 
high energies and is negligible for the bottom production. We have shown 
how to generalize the Saturation Model to the case when one of the pho- 
tons fluctuates into light vector mesons. It was found that this component 
yields a significant correction of about 30-40 % for inclusive charm produc- 
tion and 15-20 % for bottom production. We have shown that the double 
resolved component, when both photons fluctuate into light vector mesons, 
is nonnegligible only at very high energies, both for the charm and bottom 
production. 

We have shown that the production of cc pairs (the same is true for bb) is 
not completely of perturbative character and involves both short- and large- 
size contributions. The latter as nonperturbative are unavoidably subjected 
to some modeling. Present experimental statistics do not allow to extract 
cross section for the 77 — > bb reaction and therefore it is not clear where the 
observed deficit reside. It is not excluded that the apparent deficit of bottom 
quarks may reside at photon-photon energies close to threshold. This is a 
region where the underlying physics was never carefully studied. We have 
made a crude estimate of the subasymptotic quark-antiquark annihilation 
component to 77 — > bb. Although very small at high PV 77 its contribution to 
the e + e _ — > bb reaction was found to be not completely negligible. Adding 
all contributions considered in the present analysis together removes a huge 
deficit observed in earlier works on e + e~ — > bbX. 

The present analysis is based on the leading order impact factors. It 
would be desirable in the future to perform complete next-to-leading order 
calculation of heavy quark/ ant iquark production in the fc^-factorization ap- 
proach. We expect that calculating photon impact factors consistently up to 
next-to-leading order |E5] may be crucial for heavy quark production. 

Finally we have discussed a possibility to distniguish experimentally the 
different mechanisms discussed in the present paper by measuring heavy 
quark - antiquark correlations. This suggestion requires, however, further 
detailed studies of the Monte Carlo type, including experimental possibilities 
and limitations. 
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Table 1: Cross section for different contributions for charm and bottom pro- 
duction in nb and pb, respectively. 



W (GeV) 


direct 


QQ SM 


2Q2Q SM 


SR SM 


DR 


20 


1.64 


10.12 


0.79 


5.88 


0.028 


50 


0.37 


16.85 


1.35 


9.72 


0.21 


fOO 


0.11 


24.73 


1.98 


14.05 


0.67 


200 


0.034 


35.76 


2.90 


20.16 


1.78 


500 


0.0065 


58.67 


4.78 


32.06 


5.02 


20 


53.47 


301.9 


0.37 


73.06 


0.3173(-5) 


50 


14.77 


566.3 


4.16 


132.2 


0.4716(-3) 


fOO 


4.94 


840.2 


6.38 


196.0 


0.3681(-2) 


200 


1.56 


1228.0 


9.38 


287.3 


0.01719 


500 


0.32 


2047.0 


15.56 


475.9 


0.08679 



Table 2: Cross section in pb for e + e — > bbX for LEP2 averaged energy W t 
= 190 GeV. 



direct 


bb SM 


2b2b SM 


SR SM 


sum 


L3 


OPAL 


1.21 


6.1-7.4 


0.034 


1.92 


9.3-10.6 


13.1 ± 2.0 ± 2.4 


14.2 ± 2.5 ± 5.0 
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Figure 1: The cross section for 7 + p — > ccX as a function of center-of- 
mass photon-proton energy. The dotted line is the result obtained with the 
Saturation Model, the dashed line includes kinematical threshold and the solid 
line includes in addition a suppression by the factor (1 — x c ) 7 . The thick 
dashed-dotted line was obtained in the collinear approximation and the thin 
dash-dotted line represents the LO VDM contribution. 
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Figure 2: The cross section for 7 + p — > bbX as a function of the center-of- 
mass photon-proton energy. The meaning of the curves is the same as in the 
previous figure. The experimental data point is from j\L$. 
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Figure 3: Energy dependence of the main (see Eq.fi^)) SM contribution to 
17(77 — > ccX) (left panel) and 17(77 ~~ > bbX) (right panel). The dashed line 
corresponds to the prescription ^Tl\ ) and the solid line to the prescription 
(|73|j. The dotted line is the 2-gluon exchange model result with a s = 0.35. 
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Figure 4: The correction factor R c / a as a function of W for cc (solid) and bb 
(dashed) production. 
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Figure 5: Energy dependence 0/17(77 — > 2c2c) (left panel) and (7(77 — > 2626) 
(right panel). The meaning of the curves here is the same as in the previous 
figure. 
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Figure 6: R^^iW) = for charm (left panel) and bottom (right panel) 

production as a function of photon-photon energy. 
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Figure 7: ^4 map o/ 



dp\dp2 



at W = 100 GeV for the first (left-top panel) 



and the second (right-top panel) photon fluctuating into cc. For comparison 
analogous map for light quark- antiquark pairs (left-bottom panel) and for the 
case when both pairs consist of charm quarks / antiquarks (right-bottom panel) 
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Figure 8: Hadronic single-resolved contribution to 77 — > QQX for cc (left 
panel) and bb (right panel). The dashed line corresponds to the standard 
collinear calculation as described in the text, the solid lines correspond to 
the present Saturation Model calculations without (thin solid) and with (thick 
solid) the inclusion of the suppresion factor Ssr as obtained from naive count- 
ing rules. 
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Figure 9: Different contributions to the inclusive charm (left panel) and bot- 
tom (right panel) production in the Saturation Model. The long-dashed line 
represents the dipole-dipole contribution as proposed in J5|/ with the modifica- 
tions as described in the text, the dash-dotted line the single-resolved contri- 
bution calculated as described in the text and the lower dashed line the 2Q2Q 
contribution. The dotted line corresponds to the direct contribution calculated 
according to ^J\J . The double-resolved contribution is shown by the gray solid 
line. The experimental data for inclusive c/c production are from Ref.fi2Qj. 
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Figure 10: da/dW^ for e + e~ — > ccX . The main SM component is shown by 
the dashed line, the single-resolved component of SM by the gray solid line 
and the direct component by the dotted line. The thick solid line is a sum of 
all components. The experimental data are from WW . 
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Figure 11: The dependence of the integrand of on the photon-photon 

energy for xp = (solid) and xp = ± 0.5 (dashed) for the bb production 
for the direct mechanism (left panel) and in the dipole-dipole scattering in the 
Saturation Model (right panel) with the present prescription for the energy 
dependence of the dipole-dipole cross section. In this calculation W ee = 190 
GeV. 
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Figure 12: Subasymptotic quark- antiquark annihilation component (solid) 
versus the main asymptotic component of the Saturation Model (dashed) and 
the QPM component (dotted) for charm (left panel) and bottom (right panel) 
production as a function of photon-photon energy. 
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QQ correlations 




Figure 13: The expected locii in (xq,xq) space of different mechanisms con- 
sidered in the present analysis. SRi ^SR,2 means that the first/second pho- 
ton was transformed into vector mesons and DR means that each of the both 
photons was transformed into a vector meson. The dashed circle is the locus 
corresponding to the pairs emitted from the middle of the gluonic ladder (not 



